Irreversibility and self-organisation 
in hydrodynamic echo experiments. 
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We discuss the reversible-irreversible transition in low-Reynolds hydrodynamic systems driven by 
external cycling actuation. We introduce a set of models with no auto-organisation, and show that 
a sharp crossover is obtained between a Lyapunov regime in which any noise source, such as thermal 
noise, is amplified exponentially, and a diffusive regime where this no longer holds. In the latter 
regime, groups of particles are seen to move cooperatively, yet no spatial organization occurs. 
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The behavior of a system that retraces its steps af- 
ter a reversal in its dynamics - the echo protocol - has 
been of long-standing interest. From a fundamental per- 
spective, the question of how macroscopic irreversibility 
arises from the reversible dynamics of the microscopic 
components led to one of the central debates on the foun- 
dations of statistical mechanics between Loschmidt and 
Boltzmann [l[ . Loschmidt argued that reversing the ve- 
locities of all the particles in a box should allow for the 
system to return to its initial position, thus invalidat- 
ing the notion of an arrow of time. In later years, echo 
experiments became practical tools as Nuclear Magnetic 
Resonance and Neutron Spin-echo techniques to probe 
irreversible microscopic events in soft materials, and has 
also, as the name would suggest, numerous applications 
in acoustics Q. From a numerical perspective, Levesque 
and Verlet showed how rounding off errors suffice to de- 
stroy reversibility in classical Hamiltonian systems Q . A 
conceptually closely related situation, in the field of hy- 
drodynamics, is the classic Taylor experiment to denmon- 
strate the prefect reversibility of viscous fluid flows [4j. 
It consists of adding a drop of dye to a viscous fluid in 
the gap between two concentric cylinders. The drop is 
then strongly stretched by turning the inner cylinder. By 
subsequently imposing the reverse rotation, mixing is not 
enhanced: on the contrary, the initial spherical shape of 
the drop is recovered. A rather spectacular elaboration 
of this experiment has been recently performed by Pine 
et al. [f|, who substituted the ink droplet with a high 
volume fraction of non brownian beads. Although shear 
flow induces an effective diffusion of the particles in the 
suspension due to hydrodynamic interactions 0, 0] , the 
complicated particle trajectories thus generated during 
half a cycle should in principle be retraced in the second 
half. By measuring the net particle displacement after 
each back and forth cycle, Pine et al. uncovered a re- 
markable transition from a reversible situation in which 
particles do retrace their steps to a regime obtained above 
a critical strain amplitude, in which reversibility is lost. 

Two different but related questions immediately arise: 



(i) what is the origin of irreversibility, and (ii) is the 
transition a sharp one, and if so, what is its nature? 
In a recent paper, Corte et al proposed that the key 
point could be the close particle encounters, wehre non- 
hydrodynamic - and hence non-reversible - contact inter- 
actions (Van der Waals, mechanical friction) act. In one 
stroke this provides both an explanation for irreversibility 
and a mechanism for a sharp transition: because particles 
move away from situations that generate diffusion, they 
self-organise in subsequent cycles into a configuration in 
which they avoid encounters, and once this is done irre- 
versible diffusion stops. Above a certain percolation-like 
transition the system is no longer able to find such a con- 
figuration without collisions, and irreversibility persists 
in time. This scenario was further studied analytically 
by Menon and Ramsawami [^] , who made the relation to 
percolation quantitative. 

In this paper we test the opposite scenario: we 
study models that, by construction, cannot lead to self- 
organization. We assume that small irreversible pertur- 
bations, such as Brownian motion, exist all along and are 
amplified by chaoticity. This will take us to a situation 
that is very close to that of echoes in classical Hamilto- 
nian systems. We shall show that a sharp increase of par- 
ticle diffusion can be observed in such an " hydrodynamic 
echo experiment" even if no real dynamic phase transi- 
tion and no self-organization process occurs. Although 
it is likely that the contact-induced self-organization sce- 
nario is relevant for the experiments in Ref. [f| , we believe 
that other experimental situations (and perhaps even the 
numerical hydrodynamic simulation in [5j) may well cor- 
respond to simple chaotic amplification. 

Consider N force free particles immersed in a vis- 
cous incompressible fluid bounded between two concen- 
tric cylinders, the inner one being able to rotate at fre- 
quency ui(t). For spherical particles and creeping flows, 
the particle velocities are linearly related to ui [lfj : 

x a = M a (xi, ...,Xjv)w, (f) 
where the x a and x a are the particle center positions and 
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velocities. We shall henceforth set \oj\ to one. The mo- 
bility coefficients M a are complex vectorial functions, ob- 
tained in principle by solving the stationary Stokes prob- 
lem with the instantaneous boundary conditions, and 
eliminating the angular velocities using the fact that the 
particles are torque-free. The evolution is reversible, if 
we make half a period T/2 with u> and subsequently T/2 
with — ui, each particle retraces its steps over a cycle as 
in Taylor's experiment. Next, consider the effect of a 
noise i] a (t) of small amplitude, e, acting on the particle 
a that we shall assume white but in general with spatial 
correlation (r] a (t)Vb{t')) = 2e 2 S(t — t')M a b(x a , Xb). In the 
specific case of the thermal noise, M a b(x a , Xb) is equal to 
the mobility matrix, such as to respect the Fluctuation- 
Dissipation theorem ll| . The equation of motion for the 
N particles is 



Indeed, during each half-cycle, P evolves according to the 
Fokker-Plank equation , 



+M a (x 1 ,...,x JV ) +%(*), if Q<t<T/2 
-Af„(x 1 ,...,x w )+»; a (t), if T/2<t<T 



(2) 



For e ^ 0, given that the noise is different in the two 
semicycles, the trajectory does not retrace its steps ex- 
actly. How much does a trajectory deviate from the ini- 
tial position after a cycle? Importantly, the invariance 
upon time reversal of Eq. [2] tells us that this question is 
equivalent to asking how much two trajectories, starting 
from the same position at mid-cycle diverge under the ef- 
fect of different noise realisations. Let us first make this 
discussion for small noise, and for large cycle times. In 
this limit, we can express the separation, A, of the parti- 
cles between the two trajectories in the 3iV dimensional 
phase space in terms of the Lyapunov exponents [12] 



A oc e AtraJ T > 2 



(3) 



where A traj is the largest Lyapunov exponent. An impor- 
tant remark that should be made is that, because A tra j is 
a function of the trajectory, depending on the specifics of 
the driving, it may or may not coincide with the one of 
a randomly chosen trajectory after a few cycles. For in- 
stance, the system may drift away from the regions with 
high Lyapunov exponents and converge to a subset of the 
phase space corresponding to smaller diffusivity. We will 
refer to such a phenomenon as self- organization. Clearly, 
this definition of self-organization also holds for larger 
noise levels, when the Lyapunov linearization no longer 
holds. 

We now restrain ourselves to a limiting case in which 
self-organization does not happen. Let us consider the 
limit of pointwise particles, or equivalently of very dilute 
suspensions. In this limit, the particles behave like fluid 
tracers. It thus follows from the (fluid) incompressibility 
conditions that the mobility coefficients obey: 



^V a M a = 0, 



(4) 



ab 



(5) 



with V a - = d ■ /<9 Xa . This implies that the probability 
distribution P(xi, ...,xj\j) converges to the fiat measure. 



which admits a constant function as a solution by virtue 
of Eq. [U In this case, there can be no self-organisation. 
In fact, if the initial positions where chosen with flat 
probability, the probability distribution does not evolve 
and no static one-time correlation function depends on 
the cycle number. Note that there are two independent 
issues concerning external forces and interparticle poten- 
tial interactions: (i) whether they break time reversal, 
and (ii) whether they preserve or not the flat measure. 
Needless to say, in a system with divergenceless forces 
there can be no structural phase transitions, other than 
those of equilibrium hard spheres. Furthermore, the sep- 
aration of nearby trajectories corresponds to that of an 
equilibrium unbiased system, and in particular Lyapunov 
exponents are the 'typical' ones sampled in equilibrium. 
Because of the obvious analogy with the well-studied clas- 
sical and quantum mechanical problem, we shall call this 
the hydrodynamic Loschmidt-echo situation. 

To stress the analogy, we shall in this paper substi- 
tute hydrodynamic evolution, Eq. [2l by simple Hamil- 
tonian dynamics, in the presence of random noise, such 
that it conserves the microcanonical distribution. Our 
first example is very similar to the original Loschmidt ex- 
periment: we consider particles in two dimensions with 
power law r -3 interaction, perturbed by small energy- 
conserving noise (the choice of a long range interaction, 
by analogy with hydrodynamic coupling, is intended to 
limit the role of particle encounters). The second exam- 
ple, which allows us to go to large sizes and times, is a 
system of coupled simplectic maps [IH, 14 1 . 

Consider first the two-dimensional system. We per- 
form direct followed by time-reversed evolution, after a 
velocity-reversal in mid-cycle. If the system is started in 
a thermalized microcanonical configuration, then config- 
urations are statistically distributed in the same way at 
all times. Figure Q] shows the average quadratic dispace- 
ment in one cycle in terms of the cycle time, for two values 
of noise amplitude. The shape of the curve is very sim- 
ilar to the one reported for sheared suspensions and for 
superconducting vortices @, [lj| . The mean squared dis- 
placement first increases exponentially with T and then 
saturates to a constant value. 

This behavior can be explained as follows. At short 
periods, back and forth trajectories differ only slightly, 
and the Lyapunov linearization applies. Because, in the 
thermodynamic limit, there is a stable Lyapunov den- 
sity function and two subsequent Lyapunov exponents 
differ by 0(1/N) there are then O(N) largest Lyapunov 
directions which contribute to the instability, since the 
separation the y in duce are indeed comparable at any fi- 
nite time, see 12|, [l6[ . This implies that the separation 
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FIG. 1: Diffusivity D versus period T for the 2D system com- 
posed of JV = 250 particles. Over each half-cycle, the particles 
positions evolve according to: x a = — "£2 b at\xb — x a | 4 + e a , 
where the e a are uncorrelated Gaussian white noises of vari- 
ance e and a = 2 10 -6 . The box size is set to 1 and the energy 
is chosen equal to E ss 87. The diffusivity is averaged over 100 
cycles. The open symbols correspond to the zero additional 
noise limit, e = 0, yet the trajectories are not reversible. This 
is due to the (irreversible) numerical rounding off [3;]. The 
filled symbols correspond to e = (l/\/2)10 -7 . Inset: same 
plot in log scale. 



per particle is stable in the thermodynamic limit, ev- 
erything else (period, particle and energy density) being 
kept equal. At times T* such that ee AT is a sizable 
fraction of the interparticle distance, the linearization 
breaks down [l7|. Moreover, for T » T* the system 
completely loses memory in a cycle, and the separation 
becomes diffusive rather than exponential. Note, first, 
that the apparent transition time T* ~ A -1 lne depends 
very weakly on the irreversible noise amplitude, and, sec- 
ond, that it has a well-defined (large N) thermodynamic 
limit. In turn, the main features of the fluctuations of 
the stroboscopic positions with T are insensitive to the 
specific process that breaks the time reversal symmetry. 
Note also that the saturation of the actual particle self- 
diffusion occurs when the typical displacement of order 
of the box size. If we now make plots with the y-axis in 
linear scale and normalised by the saturation, the curve 
will have a sharp inflection, Figure 1. In this sense, and 
probably only in this sense, there is a sharp transition 
without self-organisation. 

In figure [5] we show the same results for a toy model: 
a time-discrete Hamiltonian system of N coupled maps 
constructed by the composition of the two following one 
step iterations and their time-reversed [l4j |: 



P'i 



: Pi 



Pi + 



)dl 




T E;=iSin [27r( g ?-g})] 



(6) 



mod 1 (7) 



FIG. 2: Diffusivity, D, versus period, T, for N = 7000 cou- 
pled maps. The diffusivity has been averaged over 500 cycles. 
Open (resp. filled) symbols correspond to a noise amplitude 
e = (1/^2)10~ 12 (resp. e = (lV2)10~ 10 ). Inset: same plot in 
log scale. 



where q[(t) = %(t + 1), p'^t) = pi(t + 1). £ Pi and £ qi are 
Gaussian white noise of variance one. Again the fluctua- 
tions of the stroboscopic position in the (q,p) plane dis- 
play the same features. Note that for this simple model 
we achieved N = 7000 and have checked that the curves 
for different (large) N coincide within numerical error 
and did not observe any sharpening of the crossover. 
Figure [3] shows the motion of particles just above the 
mobility time. Surprisingly enough, even in the absence 
of self-organization the particle system seems not entirely 
devoid of spatial features. At a given time, there are three 
types of particles: almost immobile, those that move back 
and forth along essentially one-dimensional trajectories, 
and those that have a more typical diffusive character. 
This latter group - which well above the crossover be- 
comes dominant - looks spatially correlated, suggesting 
that these are the particles that interact the most along 
trajectories. In Fig. |4]we show histograms of particle dif- 
fusivity in one cycle. The distribution seems exponential, 
and there are relatively many highly diffusive particles 
at any time. The observation of some degree of dynamic 
heterogeneity does not contradict the fact that static cor- 
relations are those of the equilibrium system, since this 
involves a purely dynamic correlation, the spatial corre- 
lations between particle displacements in a cycle. 

The scenario discussed here is by construction the op- 
posite of the one considered in Refs [1, Q , which involves 
a dramatic change in the stationary end-of-cycle distri- 
bution near the transition, completely absent here. In 
Ref. [H, the time for reaching the stationary value of 
diffusivity was plotted, showing a peak near the transi- 
tion strain, an evidence of organization. The situation 
could however be mixed: consider a system with per- 
fect time-reversal broken by a constant Brownian noise, 
so that the flat distribution is preserved. Suppose we 
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FIG. 3: Stroboscopic plot of the particles position after each 
cycle for the 2D system of interacting particles. Plot after 
100 cycles for a period T = 2, same parameters as in Figure 1 
(open symbols). 




-0.05 



FIG. 4: Histograms of particles diffusivity for T = 2 (dark red 
circles), T = 2.2 (light red circles) and T = 2.4 (open circles). 
Same parameters as in Figure 1 (open symbols). Inset: Same 
histograms for the coupled maps with e = (l/\/2)10~ 12 and 
T — 28 (dark blue squares), T — 30 (light blue squares) and 
T = 32 (open squares). 



switch on forces that violate the uniform distribution: 
this indeed is what happens in the hydrodynamic case as 
the volume fraction is increased. The system will now 
self-organize into a distribution that, unlike the case we 
studied here, will depend on strain, noise level and par- 
ticles volume fraction. If the effect is weak, we do not 
expect that a sharp transition will immediately arise, 
so that from this point of view the scenario is still the 
one we discussed here, in spite of having a certain de- 
gree of auto-organisation. In order to further empha- 
sise the experimental relevance of the hydrodynamic echo 
scenario, we compute the crossover shear amplitude 7* 
which would correspond to the experimental conditions 
reported in [jj Q. The noise source is supposed to be 
the weak Brownian motion of the suspended particles, 
e is thus approximated by the distance over which a 
single particle diffuses over half a cycle: e = V DT, 
D = k-oT/ (Qirrja) with 77 = 3Pa.s the fluid viscosity and 
a = 115/im the particles mean radius [s] (Note that we 
have neglected hydrodynamic correlations of the noise). 
The interparticle distance is of order of A ~ a/cj) 1 ^ 3 . To 
compute 7*, we also need the Lyapunov exponent of the 
concentrated particle suspension. This value has been 
computed numerically in [l8|], for <j> — 0.4: A ~ O.67 
where 7 = 2n~f/T for sinusoidal cyclic shear. We thus 
infer a critical shear amplitude of: 7* ~ 3 which is fairly 
close to the experimental value: 7* xp ~ 1 given the crude- 
ness of of our numerical estimates. This predictions sug- 
gest that there are very probably many situations in- 
termediate between the pure Loschmidt echo crossover 
discussed in this paper and percolation-like transition. 
In such cases, the relative importance of the two effects 
cannot be solely revealed by the measurement of the ef- 
fective diffusivity of the particles and may be harder to 



assess experimentally. 

J.-L. Barrat, P. Chaikin, L. Cipelletti, L. Corte, G. 
Menon, D. Pine, S. Ramaswamy and especially D. Levine 
are acknowledged for valuable comments and discussions. 
G. D. was supported by a CONIC YT grant. 



[5 
[6 
[T 
[§ 
[9 
[10 

[11 

[12 
[13 
[14 

[15 



R. H. Swendsen Am. J. Phys. 76, 643 (2008) 
M Fink, et. al., Rep. Prog. Phys. 63, 1933 (2000) 
D. Lesveque and L. Verlet, J. Stat. Phys. 72, 519 (1994) 
G. I. Taylor, National Committee for Fluid Mechan- 
ics Films, Education Development Center, Newton, MA 
(1966) 

D. J. Pine, J. P. Gollub, J. F. Brady, A. M. Leshansky, 
Nature 438, 997 (2005) 

E. C. Eckstein, D. G. Bayley and A. H. Shapiro, J. Fluid. 
Mech. 79, 191 (1977) 

D. Leighton and D. Acrivos J. Fluid. Mech. 177, 108 
(1987). 

Laurent Corte, P. M. Chaikin, J. P. Gollub and D. J. 
Pine, Nature Physics 4, 420 (2008) 

G. I. Menon and S. Ramaswamy, cond-mat 
larXiv:0801.388TV l (2008) 

J. Happel and H. Brenner, Low Reynolds Number Hydro- 
dynamics, Springer (1981) 

P Pusey, in Liquids, Freezing and the glass transition 
Edited by JP Hansen, D Levesque and J Zinn- Justin, 
Elsevier Amsterdam (1991) 

G. Veble and T. Prosen, Phys. Rev. E, 72, 025202(R) 
(2005) 

T Konishi and K Kaneko, J. Phys. A: Math. Gen. 25, 
6283 (1992) 

E. G. Altmann and H. Kantz, Europhys. Lett. 78, 10008 
(2007) 

N. Mangan, C. Reichhardt, and C. J. Olson Reichhardt, 



■5 



Phys. Rev. Lett. 100, 187002 (2008) Rep. 356, 367 (2002) 

[16] R Livi, A Politi and S Ruffo, J. Phys. A: Math. Gen. 19, [18] G. Drazer, J. Koplik, B. Khusid, A. Acrivos, J. Fluid. 

2033 (1986) Mech. 460, 307 (2002) 

[17] G. Boffetta, M Cencini, M Falcioni and A Vulpiani Phys. 



